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Abstract. Since Fermions are based on anti-commutation relations, their entanglement 
can not be studied in the usual way, such that the available theory has to be modified 
appropriately. Recent publications consider in particular the structure of separable and of 
maximally entangled states. In this talk we want to discuss local operations and entangle- 
ment distillation from bipartite, Fermionic systems. To this end we apply an algebraic point 
of view where algebras of local observables, rather than tensor product Hilbert spaces play 
the central role. We apply our scheme in particular to Fermionic Gaussian states where the 
whole discussion can be reduced to properties of the covariance matrix. Finally the results 
are demonstrated with free Fermions on an infinite, one-dimensional lattice. 



1. Introduction 

Entanglement distillation is one of the most fundamental processes of 
quantum information processing [I]. It is an integral part of many protocols 
and devices like quantum repeaters and it provides important procedures to 
measure the entanglement content of a given physical system. In the usual 
setup it is assumed that two distant parties - Alice and Bob - share a large 
number of (weakly) entangled pairs of particles, and the task is to generate 
a (possibly small) amount of highly (or even maximally) entangled pairs by 
means of local operations and classical communication. In this context it is 
implicitly assumed that the particles shared by Alice and Bob are distinguish- 
able: Firstly there is a clear distinction between particles controlled by Alice 
and those controlled by Bob. Secondly many distillation protocols require 
the selection of particular pairs by Alice and Bob (e.g. to perform filtering 
operations). Hence even the local particles needs to be distinguishable. 

At a first glance it therefore appears to be completely impossible to talk 
about entanglement distillation with fermions, because the latter are by def- 
inition undistinguishable. A second thought reveals, however, very quickly 
that this is only an apparent difficulty. The only thing we have to drop is 
the focus on particles. Instead we have to consider setups where Alice and 
Bob control independent subsystems (usually distinguished by their position 
in space) of a larger physical system. A typical example is a Fermi gas from 
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which Alice and Bob want to distill entanglement by using only operations 
(together with classical communication) which are localized in spatially sep- 
arated regions (e.g. Alice's and Bob's laboratories). 

Mathematically such a situation is most easily described in terms of op- 
erator algebras. In other words, instead of using tensor products of Hilbert 
spaces, we describe bipartite systems by specifying which local observables 
are measurable by Alice and which by Bob. This approach is successfully 
applied to the study of entanglement of infinite degrees of freedom systems 
[111 [7] and for the analysis of separable [SI [5] and maximally entangled |9j 
states of Fermionic systems (cf. also the references in [5] for more literature 
on Fermionic entanglement). 

The purpose of this paper is to study entanglement distillation in the same 
framework. In this context we will show that distillation from Fermions can 
be treated basically in the same way as ordinary distillation with only some 
small changes which mainly arise from the emergence of super selection rules. 
In addition we will present an explicit scheme which can be applied to any 
quasi-free state and which allows explicit calculations (e.g. of distillation 
rates) for fairly large systems, such as (subsystems of) infinite quasi- free 
lattice models. 



Let us start with a short look on standard distillation techniques. Hence 
assume that Alice and Bob share N d- level systems in the joint state p® , 
where p denotes a density matrix on the Hilbert space TL <8> Ti, 7i = C d . 
To generate maximally entangled qubit pairs from these resources they can 
proceed as follows: 

1. Look for (and drop) unentangled subsystems. Mathematically this 
means to find a unitary U : Ji — > H\ (8> K., H\ = C^ 1 , and to apply 
the transformation T\(p) = tr^(UpU*). The final state p\ should con- 
tain (almost) as much entanglement as the original p. In other words: 
JC and U has to be chosen appropriately. 

2. Find a maximally entangled state ip G H\ £3 H\ such that 



The best choice would be to take the tp which maximizes this fidelity. 
To be successful here the first step is in many cases mandatory, be- 
cause a state can be entangled without satisfying inequality dU) for any 
maximally entangled tp. This can be easily seen if we choose d = 2d 



2. Entanglement distillation 
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and 

p=|*WI®4 <t> = 2-^(100) + |11» (2) 

since we get sup^{ip, pip) = 1/d 2 , which can be arbitrarily small (if d is 
big enough) although p is distillable. 

3. Consider the group 

G$ = {U A ® U B I U A , U b e U(di), C/a 8> W = */>}, (3) 
and average over it. This leads to the twirl operation given by 

T 2 (pi) = I UpU*dU, (4) 



where dU denotes the Haar measure on Gv,. 
4. The output of the channel T 2 is an isotropic state 

.0 

with i? G [— (df — I)" 1 ! 1] given in terms of the fidelity / = (ip,pip) by 



P2 = T 2 p 1 = tf|V>M + (1-^)12 (5) 



From Equation ([T]) it follows immediately that p 2 again satisfies the 
inequality 

(^p 2 ^)>d^ 1 , (7) 
and therefore it is distillable [6]. 

5. Now we can continue with standard techniques for isotropic states 
which provide us with a number M of (almost) maximally entangled 
qubit pairs; cf [12] and the references therein. 



3. Bipartite Fermionic systems 

To study entanglement of Fermionic systems the usual framework which 
relies on tensor product Hilbert spaces is too narrow because indistinguisha- 
bility and anti-commutation relations has to be taken into account. Instead, 
it is more appropriate to describe the splitting of the overall system into two 
subsystems in terms of observables algebras. This approach was successfully 
applied in particular to infinite degrees of freedom systems [111 [7] . For our 
purposes a simplified (finite dimensional) approach is sufficient. 
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DEFINITION 1. A bipartite quantum system consists of a Hilbert space and 
two C*-algebras A, B C B(7i) which commute elementwise (i.e. [A, B] = 
Gi,VBG B). 

Selfadjoint elements of A and B describe the (projection valued) observ- 
ables of the given system which can be locally measured by Alice and Bob 
respectively. The default setup can be recovered if we have 

A = B(H a )®*b, B = 1 a ®B{H b ), M=B{H a ®H b ) (8) 

in terms of two Hilbert spaces Ti.A,T~(-B- Note, however, that we have neither 
assumed that A and B together generate B(H) nor that B is the commutant 
of A (in contrast to [7])- Therefore beside ([8]) other realizations of bipartite 
systems are possible, even if H is finite dimensional. A particular example 
arises, if TLa/b decomposes into a direct sum TLa/b = ^a/b ® ^A/B' an< ^ ^ 
we define 

A = {B(H\) © B{H A )) <g> 1 B , B = 1 B ® (B{H\) © B(H^)) . (9) 

This is - as we will see - exactly the situation we have to study for a system 
consisting of a finite number of Fermions. 

To explain the latter remark consider the Hilbert space /C and the corre- 
sponding antisymmetric Fockspace H = J--(JC). For each h, f E K we can 
define the usual creation and annihilation operators c*(h) and c(/) on Ti., 
which satisfy the canonical anti-commutation relations ({ • , • } denotes the 
anti-commutator) 

{c(f),c(h)} = {c*(f),c*(h)} = 0, {c*(h),c(f)} = (h,f)l. (10) 

In some cases it is more convenient to combine c and c* in one operator (this 
is called the self-dual formalism [2j [3] ) : 

B(f,h) = c(f) + c*(h), T(f,h) = (KJ), (11) 

where (~) denotes complex conjugation in an appropriately chosen (and 
fixed!) basis. Now we get from (jlOf) 

{B(F 1 ),B(F 2 )} = (rF 1 ,F 2 ), B(F)* = B(TF),VF, Ft, F 2 e /C © £. (12) 

The set of operators {c(/) | / G IC} C B{H ) generates a C*-algebra CAR(/C) C 
B(7i) which is called the algebra of canonical ani-commutation relations. It 
can be regarded as the closure (in operator norm) of the algebra of polyno- 
mials in the c(/) and c*(h). 
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Now consider the parity operator on H which is given in terms of the 
number operator N by 6 = (—1) N . It acts as the identity on the subspace 
7i + C H of vectors with an even particle number and as minus the identity 
on the complementary subspace TL~ . We write P^ for the corresponding 
projections and get 

P ± :H^H ± , 9 = P+-P-. (13) 

The elements of CAR(/C) which commute with 9 are called even elements 
and they form the even subalgebra 

CAR+(/C) = {A G CAR(/C) | [A, 9} = 0} (14) 

of CAR(/C). It can be regarded as the (closure of) the algebra of even poly- 
nomials in c(/) and c*(h). 

If IC is finite dimensional CAR(/C) coincides with B(7i), i.e. it is a full 
matrix algebra. The even subalgebra however is given by 

CAR+(/C) = B(H + ) © B(H~). (15) 

This can be easily seen from the fact that a product of an even number of 
creation and annihilation operators can change the particle number only by 
a factor of 2. 

The next step is to decompose IC into an "Alice" and a "Bob" subspace, 
i.e. IC = JCa © ICb- Then we can associate to K>a/b the corresponding 
Fockspaces Ha/b = ^-{^a/b)i an d also the CAR-algebras CAR,(1Ca/b) an d 
CAR+(X^4/b). Obviously we have 

H = T-{K A ) ® F-(K.b) =H A ®n B , (16) 

and similarly CAR(/C) is isomorphic to CAR(/Ca) ®CAR(/Cs) if we consider 
the spatial tensor product. The corresponding isomorphism satisfies 

c(fA(B0)^c A (fA)®K B , c(O0/ B ) ^9 B ®c B {F B ), (17) 

where the ca/b(J 'a/b) denote the annihilation operators and 9a/b the parity 
operators on TLa/b- The latter are related to the global parity 9 by 

9 = 9 A ®9 B . (18) 

Equation fjlT[> shows that CAR(JCa/ B ) can not both be embedded as tensor 
factors into CAR(/C) without violating the anti-commutation relations. The 
construction given in f)17j) is therefore called the twisted tensor product of 
CAR(/C A ) and CAR(/C B ). 



M. Keyl, Distilling entanglement from Fermions 



6 



The discussion of the last paragraph has shown that two operators A E 
CAR(ICa), B E CAR(/C#) do not commute in general and therefore these 
two algebras can not be chosen as the observable algebras A and B. If we 
choose, however, A and B to be even elements we immediately get from (|17p 
that [A, B]=0 holds. Hence 

,4 = CAR + (/C A ), £ = CAR + (/C b ) (19) 

is an appropriate choice for the local observable algebras of Alice and Bob. 
Together they define a bipartite Fermionic system in the sense of Definition 
f3TI If JCa and K, b are finite dimensional we can insert Equation (115f) for A and 
B and we recover the example already given in ([8]). Note in addition that 
A C B{T{a) and B C B(Hb) holds, while the full CAR algebras fail to have 
this property - by virtue of Equation (|17|) . Hence it is reasonable to consider 
~Ha (Hb) as Alice's (Bob's) Hilbert space. 

4. Local operations 

On top of the scheme described in the last section all basic notions of 
entanglement theory can be reconstructed. This was done for separable states 
in [U [5] and for maximally entangled states in [9]. To discuss entanglement 
distillation we need the concept of a local operation which can be defined as 
follows (cf. 0): 

DEFINITION 2. Consider two bipartite systems Aj,Bj C B(Hj), j = 1,2. 
An operation (i.e. a completely positive map) T : B(7ix) — * B^H.2) is called 
local if 

T{A X ) c A 2 , T{Bx) c B 2 (20) 

and 

T{AB) = T{A)T{B) V-A E Ax, V-B E B\ (21) 

holds. 

In the standard framework ([8|) with finite dimensional Hilbert spaces 
T~(-A,7~tB (or if we assume that the operation T is normal) this definition 
coincides with the usual one. Note that the factorization condition (|21l) is 
needed to make this statement true [?]■ 

To generalize the distillation protocol from Section [2] only a few special 
local operations are needed. They are summarized in the following list. 

• Local unitaries. The easiest case is a local unitary transformation A i— > 
U*AU with U = Ua®U b and Ua (Ub) unitary on Ha (He)- It is easy 
to see that U* A AU a = A is equivalent to U aH-a = %A or ^ A^-a = ^A- 
A similar statement holds for Ub- 
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• Local Bogolubov transformations. A special case of the previous exam- 
ple arises if Ua and Ub are related to unitaries ua, ub on JCa © K-A-, 
JCb © ICb by 

U*B(F)U = B(uF), with U = U A ®U B , u = u A ® u B , (22) 

where B(F) is the operator introduced in Equation (jlip . It is easy to 
see that (i22~j) is only possible if 

Pur = u (23) 

holds (cf. again (jlip ). The condition (|23(1 is on the other hand sufficient 
for the existence of a unitary U satisfying (|22p for a given u. 

• Discarding subsystems. Consider now decompositions K-a/b = ^U/B,l© 
K-a/b,2 °f K-A and /C#. If we denote the Fockspaces of Kj = K.A,j®K>B,ji 
j = 1,2 by Ttj we get a decomposition of W into a tensor product 
H = 7ii © H2. If we perform the partial trace over say H2 we get 
the Fockspace Tii and the corresponding reduced observable algebras 
Ai,Bi C BiTii). In this way the partial trace becomes (the Schrodinger 
picture version of) a local operation between two bipartite Fermionic 
systems, which discards the modes belonging to K,2. 

• A joint parity measurement is described by the PVM 

P jk = P j A®P k B , j,k = +,- (24) 

and the corresponding von Neumann-Liiders instrument {Pa/b d en °t e 
the projections to the even/odd subspaces of Ha and Tis', cf. Section 
[37]) . For a system in the state p the probability to get the outcome j, k 
is pi k and the corresponding output state is pP k : 

pjk „pjk 

p* = tr(pi k p), pi k = P r . (25) 

The projections P^ k commute with all A G A and all B £ B. Therefore 
parity measurements can be done without disturbing the system. This 
implies immediately that the state p can not be distinguished from the 
mixture Yljk'P' r Within a distillation scheme this instrument can 
be used to perform local filtering operations; e.g. Alice and Bob can 
decide to drop the whole system if their local parities are different and 
to keep it otherwise. If we set 

JC A = JC B = C d ^Ti\ = Ti^ = Ti^ = Ti B = C D , D = 2 d -\ (26) 
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and ignore the value of the parities (apart from j = k) we get a (non- 
unital) local operation which transforms a bipartite Fermionic system 
into a pair of D— level systems in the state 



p++ p++ _|_ p p 
p++ _|_ p 

5. Distilling from Fermions 



(27) 



Let us now adopt the general distillation scheme sketched in Section [2] 
to the Fermionic case. To this end we will use throughout this section the 
assumptions made in Equation (|26p . which implies in particular that Q 
holds. In addition, consider two maximally entangled vectors (p+ £ 7i ++ , 
tp- £ 7i and 

t/> ± = ^=(tp+±<pJ). (28) 

For each A £ A <S> B we have 

tr(A|V±>(V±|) = \ (tr(A| V+ )( V+ |) + tr(A|^_)(^_|)) . (29) 

Using the terminology from Equations (|45p and ()57[) this can be rewritten 
as: 

P ++ = P~~ = x, P + ~ =P~ + = 0, p ++ = \<p+)(<p+\, p~ = \ip-)(<p-\. (30) 

Hence Alice and Bob can not distinguish the vector states |V'±)(V'±I f rom 
themselves and from the mixture of \(p + )((p + \ with The latter is 

according to [9] a Fermionic maximally entangled state (implying in partic- 
ular that EOF is maximal). 

The only step from the list in Section [2] we have to change is the twirling, 
because averaging over the group (or G^_ ) breaks the superselection rule 
and is therefore not an allowed local operation. Instead, we have to look at 
the subgroup 

H f+ = {U A ®U B e GV + | UaAUa = A, U B BU B = B}. (31) 
The structure of this group is given by the following Proposition 
PROPOSITION 3. The group is generated by the subgroup 

H<p + , = {U A ®U B £G t p + \ [U A , 0a] = 0, [U B ,e B ] = 0} (32) 
and V = Va <8> V B given by 

V Ae%j = e~ A)j , V B e£j = e Bj (33) 
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decomposition of ' <p±, i.e. 



1 ° 1 ° 

^+ = ^E e Xi ® e B ,j » <p- = ~m E e A,j ® e B,i • ( 34 ) 



Proof. Obviously H^ +i q C -£^ + and Va (g) Vb G To show the other 

inclusion recall from the discussion of local unitaries in the last section that 
U A AU B = A is equivalent to V r A U\ = U\ (i.e. [U A ,0 A ] = 0) or U A H^ = 
TC^, and that a similar statement holds for Ub- The assumption U A (&U sip+ = 
Tp + implies in addition that [E/a>0a] = 44> [Ub,0b] = holds. Hence 
U £ is either in H^.q or it can be written as U = UV with a U £ H^.q, 
which concludes the proof. □ 

Averaging over the group leads to states which are invariant. 
Their structure is given by the following proposition. 

PROPOSITION 4. Each -invariant state a can be written as 

a = A+IV+X^+I + A_|^_)(^_| + MP + + + P~) + V-(P + ~ + P~ + ) (35) 
with 

p ++ =p— = A+ ^ A ~ p- + =p + - =fi-D 2 , (36) 

{ip± , o-ip± ) = X± + n + . (37) 

Proof. We have to determine the commutant of H^ + . To this end 

note first that H^ +t o C implies H'^ C H'^ Q . Hence consider the latter 
commutant first. By definition we have for each unitary U on TC A ®7is 

U e H^o O [U, |^ + X^ + |] = 0, [U,6 A ® 1] = 0, [U, 1®0 B ] = 0, (38) 

where we have used the fact that the factorization U = U A (g> U B is a conse- 
quence of [U,ip+] = 0; cf. [13j. Therefore H'^ + is the von Neumann algebra 
generated by \ifj + )(ip + \, 6 A (g) I and 1 (g> 0b, i.e. 

^+,0 = {\iI>+)W+\,0a ® I, I® 0b}". (39) 
By calculating all possible products of the generators this leads to 

H ip+,o = span{|'0 + )(^+|, \ip-){ip-\, |^+)(^_|, |^_)(^+|, 

P ++ ,P—,P + -,P- + }. (40) 
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The group is generated by Hip + ,o and V = Va®Vb', cf. Proposition 
ED Hence A G 

o i s i n ^ ^ commutes V. Since Va and Vb just 
exchanges the even with the odd subspace we easily conclude that 

H'^ + = S p & n{\^ + )(i; + \,\^)(^\,P ++ + P— ,P+- + P-+} (41) 



holds, which implies equation f)35[) . Equations (|36p and (|37p follow immedi- 
ately from the definition of the p^ k in (145j) and from taking traces. □ 



If we decompose the -invariant state a according to Equation (Jo 
we get 

Hence if Alice and Bob perform 9a, Ob measurements - which they can do 
without disturbing the systems - they get either with probability 

p = p ++ +p— (43) 

one of the (basically equivalent) isotropic states c ++ or a , or they get with 
probability 1—p the totally chaotic state or a In case they get 
it is distillable iff 

{<P±,o ±± <p±)>^ (44) 



holds. A straightforward calculation using Equations (|36l) . (|37|) . ([42]) and 
(|44p leads to the following proposition 

PROPOSITION 5. Consider a state H^, invariant state a. The fidelity 
f = {<P±i f ±± V 9 ±) o/ isotropic state o" ±± /rom Equation J^I| ) is given by 

p 

Hence <t ±=i= is distillable iff 

(1>+,<nl>+) + (il>-,<"l>-)>% (46) 

Let us start with a general state p and twirl over H^ + , i.e. 

a = I UpU*dU. (47) 
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Then a is invariant and Equation (|46f) is equivalent to 

(V + ,p^+) + (^-,pV-}>§ (48) 

To get a general distillation protocol for Fermionic systems we can therefore 
modify the distillation presented in Section [2] as follows: 

1. Start with N distinguishable copies of the same Fermionic system, each 
prepared in the same state (e.g. N metallic wires containing an electron 
gas). 

2. Drop unentangled subsystems. This leads to N bipartite Fermionic 
systems in the joint state p® N . Here p is a density operator on Ha®'Hb, 
which should be interpreted, however, as a state of the algebra A <g> B. 

3. Find maximally entangled states ip± as in Equation (j28j) such that (1481) 
holds. 

4. Average over the group H^ + . This leads to the -invariant state a. 

5. Make 9a, Ob measurement. If the outcome is ++ or (which happens 

with probability p) this leads to the isotropic state <r ++ or a . It can 
be treated with standard distillation techniques. 

6. Otherwise (H — , — h) we get a chaotic state which is useless for distil- 
lation. 

To find a maximally entangled state such that (j48|) is satisfied usually 
requires an optimization over all possible In general this is very difficult. 
In the next section, however, we will discuss a special class of states where 
this problem is more feasible and which provide at the same time a systematic 
way of dropping unentangled modes. 

6. Quasifree states 

Let us apply the general scheme developed in the last section to quasifree 
states. Recall that a density matrix p describes a quasifree state of the CAR 
algebra CAR(/C) if there is a bounded operator S £ B(JC) such that 

tr(pB(f l )---B(f 2n+1 )) = (49) 

n 

tr(pS(/i) • • • B{f 2n )) = ]T sign(p) n< r / P (2;-i), Sf p{2j) ), (50) 

3=1 

holds for all n G N and /& G K © K. The sum in (|50p is taken over all 
permutations p satisfying 



p(l) <p(3) < ... <p(2n-l), p(2j-l)<p(2j) (51) 
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and sign(p) is the signature of p. The covariance operator S is selfadjoint 
and satisfies 

rsr = i-s, o < s < i. (52) 

We can express the right hand side of Equation (i50"j) as the Pfaffian Pf(S') 
of the antisymmetric matrix S with matrix elements Sjk = (fj, Sfk) k > j. 
Note also that this definition works as well in infinite dimensions, although we 
will restrict our discussion in this chapter again to the case K-a = /C# = C d , 
K, = ICa © K-B- In K-a/b © K-A/B we wu l use the bases (k = 1, . . . , 2d) 



2- 1 ' 2 (\k) + \k + d)) iik<d 



(*) _ J A/B,l ~ v,-, ■ ,- ■ -- ■• - - 

A/B | (k-d) „_1/9.-/l, n l,n -r , . , V 0, 3J 



A/B,2 



2- 1 /2i(|A; - d) - |ife)) ifA;>d 



where |/c), A; = 1, ...2(2 denotes the canonical basis in C 2c( . If the decom- 
position into Alice- and Bob-subsystems is not important we can also use 
relabeled version 

e «_R } for k < 2d 

6 ~\e% d ~ k) for k > 2d (54) 

where k ranges now from 1 to Ad. The advantage of this basis is its T 
invariance. We can therefore write 

S kj = (e^ k \Se^) = (re( k \Se^), (55) 

with k,j = 1, . . . ,4<i. In the following we will identify with slight abuse of 
notation the operator S with the matrix (Sjk)j t k an d write 

S=(" AA S v AB ) = 1 -( 1a+ Jt X * Y :)- (56) 
V Sba S BB J 2 V ~tY T 1b + %Z J v ! 

These expressions should be interpreted as block matrices with respect to the 
Alice/Bob split, e.g. Sab contains all matrix elements of the form (e^ , Se^ ), 
etc. Using Equation (|52p it is easy to see that X, Y, Z are real 2d x 2d matrices, 
and that Y, Z are antisymmetric. 

For quasifree states the expressions and constructions from the last two 
sections can be given quite explicitly in terms of covariance matrices. The 
following list summarizes the most important examples (cf. |1U] for more 
details, in particular for proofs) 

• The probability p to get equal parities during a joint 0a,0b measurement 
(cf. Equation (03]) is given by 

l + (-4) d Pf(5-I/2) 

P= n ' ( 57 ) 
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A quasifree state with covariance matrix P is maximally entangled, iff 
P has (in the basis from Equation (|53p ) the form 

with a real orthogonal matrix R. The quasifree state thus given can be 
represented by a state vector ipp E TL with 

^P = ^(<^+ + ^-) (59) 
with maximally entangled vectors <p± € ® Wg. In other words ipp 



is always of the form ip + assumed in Equation (|28p . 

• The fidelity between a quasifree state ps and a maximally entangled 
quasifree state ipp is given by 

(Vp,pV>p> = Pf (2-5- P). (60) 

• The quasifree state ps can be transformed by a local Bogolubov trans- 
formation u into a normal form p$ (which is again quasifree) such 
that the off diagonal blocks —iSAB and iSpA of the covariance matrix 
S = uSu* become diagonal with positive eigenvalues. To see this con- 
sider the singular value decomposition Y = ua^yu*b °f ^ an d choose 
u = ua © up- 

The general distillation scheme described in the last section comprises 
the search for a ip+ such that Equation (|48j) holds. Since we can always 
choose ip + = ipp for some P satisfying ([58]) a good strategy is optimize the 
expression in (I60p over all such P. The following theorem treats an important 
special case (cf. [TU] for a proof). 

THEOREM 6. Consider a quasifree state ps with covariance matrix S from 
Equation $56)) . Assume that X = or Z = holds, and that Y is diagonal 
with non-negative eigenvalues (the latter can be done without loss of general- 
ity). The maximal fidelity of ps with a quasifree, maximally entangled state 
ipp arises if the basis projection P is given by 

and its value is 



2 V -tl Is 



<^,w>p> = n Hr^ (62) 

3=1 V J 

where Xj, j = 1, . . . , n denote the eigenvalues values of Z and rrij the corre- 
sponding multiplicities. 
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If the condition X = or Z = is not satisfied the optimality statement 
is in general not true. For states, however, which are already close to a 
maximally entangled, quasi free state X and Z have to be at least small 
(otherwise the condition < S < 1 is not satisfied). Hence in this case the 
choice tpp with P from (I61j) should be close to the optimum (provided Y is 
diagonalized) . Therefore the following specialization of the procedure from 
the last section should provide a reasonably good scheme for distillation from 
quasi free, Fermionic states. 

1. Consider N bipartite Fermionic systems, each of which in the same 
quasi free state p$. 

2. Choose bases for K,a®^a and Kb® Kb such that the block-offdiagonal 
part Y of S becomes diagonal (and with real positive entries). As 
already pointed out above this can be done locally by Alice and Bob 
without any communication (if S is known to them). 

3. Drop all modes except those belonging to the n highest singular values 
of Y . The number n must be chosen such that Equation (148j) holds with 
ip + = if) P where the basis projection P is (in basis which diagonalizes 
Y) of the form (I61j) . The maximally entangled state ip- is then quasi 
free as well, i.e. ip- = ipQ with basis projection 



(please check yourself). Hence the two fidelities {ip±ps>' l P±} m Equa- 
tion (148 h can be calculated with (I61j) . If the probability p is unknown 
Equation (|48p should be used with the conservative choice p = 1. 

4. Average (twirl) over the group H^ + , make a 9 a, Ob measurement and 
proceed as described in Section [571 

Let us demonstrate this scheme with free fermions (without spin) hopping 
on a one-dimensional regular lattice Z (lets call it a "wire"). They can 
be described by the CAR algebra CAR(1 2 (Z)) and the dynamics is given 
formall}0 by the Hamiltonian 



1 H is not well defined as an element of CAR(1 2 (Z)), because the sum does not converge 
in norm. It gives rise, however, to a well define derivation and therefore the notion of 
ground state is well defined too. 




(63) 




(64) 
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It admits a unique quasifree ground state ipo with covariance operator S given 
by 

where 

oo 

1 2 (Z) <g> C 2 3 F h+ e L 2 ^ 1 ) <8 C 2 , .F(F)(x) = ^ e inx F n (66) 

j=-oo 

is the Fourier transform and E £ S(L 2 (S' 1 )) the projection to the upper 
half-circle [2JS]. 




Fig. 1: Distillation rate i?(<i) /d for free Fermions on a one-dimensional lattice, 
if only two adjacent regions of length d are accessed. 

Now assume that Alice and Bob can control only two blocks of the form 
Aa = [0,d) and A# = [d, 2d) (or any joint spatial translate of them). The 
restriction to the corresponding subsystem leads exactly to the bipartite 
Fermionic system just studied. The reduced density matrix p\ arising from 
the ground state <po is quasifree and its covariance matrix in the basis (|54|) 
can be easily derived from ([65]) . 

Now we can apply the distillation protocol given above. If we choose 
to keep in step [3] only the four highest singular values of Y we get at the 
end with probability p from (I57p a qubit pair in an isotropic state a = 
f\<p)(<p\ + (1 - /)(! - \4>}(4>\)/3, with fidelity / from Equation (05]). If a large 
number of systems is available (where system refers here to a whole wire not 
to a single Fermion) and if the fidelity / is big enough we can use the Hashing 
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protocol to distill maximally entangled qubit pairs. The distillation rate, i.e. 
the number of maximally entangled pairs we get asymptotically per wire is 

m 

R = p(l-S(a)) =p(l + /log 2 (/) + (l-/)log 2 (l-/)-(l-/)log 2 (3)). (67) 

Maybe more interesting is the rate R/d of pairs we get per lattice site used. 
The result is plotted in Figured! The small zigzag noise on the graph arises 
from a slightly different behavior of the protocol for even and odd values 
for d. This is an indication that the scheme is indeed not optimal if the 
assumptions from Theorem [671 (i.e. X = or Z = 0) are not satisfied. 
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